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We report results on semileptonic B irlV decay form factors near ^max using NRQCD heavy quark and clover 
light quark actions and currents improved through 0{aa). An inconsistency with the soft pion relation /"(^max) ~ 
/fl //tt found in a previous work is confirmed, and a possible solution with nonperturbative renormalization is 
discussed. We find that f^iq^) is well described by the B* pole near Qmaxi a^nd its 1/Mb scaling is also consistent 
with the prediction of the pole dominance model. 



1. Introduction 

The B ttW form factors are relevant for 
the extraction of the CKM matrix element \Vub\ 
through the exclusive decay. While lattice QCD 
computation can cover only the region near — 
9max with reasonable statistical and discretization 
errors, it is still useful once the experiments reach 
sufficiently high statistics to measure the partial 
decay rate in the same region. In this article 
we report preliminary results of our study of the 
B ttIv form factors using NRQCD for heavy 
quark. 

2. Simulation 

We employ 0{l/M) NRQCD for heavy quark. 
For the light quark the SW clover quark action 
is used, with the clover coefficient Csw determined 
by mean field improved perturbation theory at 
one-loop order. The heavy-light vector current 
involved in the matrix element is renormalized to 
0{asa) using the one-loop calculation of Morn- 
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ingstar and Shigemitsu 0], and of Ishikawa ||] 
including mixings with higher dimensional oper- 
ators. 

The calculation of /"(gmax) was performed at 
/3=5.7 and 5.9 on 12^ x 32 and 16^ x 40 lattices 
with 232 and 222 gauge configurations, respec- 
tively. We took five different heavy quark masses 
covering m;, and four different light quark masses 
ranging from to TOs/2. For f~^{q^) we only an- 
alyzed the (3 = 5.7 data, so far. The momentum 
combinations p_B = (0, 0, 0), fc7r=(0,0,0), (1,0,0), 
(1,1,0), and PB = (1,0,0), fc^=(0,0,0) in units of 
27r/(12a) are considered at k; = 0.1369, which is 
around nis- 

3. Results for fiql^J 

Let us first present the results for /"(^max)' 
for which the soft pion theorem /°(gmax) = 
fs I fn in the chiral limit provides an impor- 
tant check of the lattice calculation. The 
1/Mb dependence of /Ms /°(9max) (a^A^s) 
/as(MB''^'))^/^° and a comparison with y/M^ 
Ib/U («.(Ms) /asiM^J'n?/^" is shown in 
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Figure 1. Comparison of ^ M b {q^a^^)ri with 
\J Mb^b I j-K"r\-, where r\ in defined as (Q!s(Afs) 
las{Mf^^)fll^°. with heaviest and 

with Ughtest fight quark mass is plotted at (3=5.7 
and 5.9. Data for ]bI J-k are taken from our pre- 
vious works: heavy clover 0] and NRQCD H. 



Fig. ||. The data for /"((/max) a-t two (i values 
show nice scaling, while a clear disagreement with 
/b/Ztt is observed, confirming the point made by 
the Hiroshima Group [||. 

The chiral (soft pion) limit is not taken for the 
points given in the plot, which is a possible reason 
of the violation of the soft pion relation. We find, 
however, that the light quark mass dependence 
of /°(9max) is consistent with a constant within 
statistical error, and a polynomial chiral extrap- 
olation in vfiq and gives a consistent result. 
With the present statistics, we are not able to fit 
our data with both ^rtiq and mg in the fitting 
function. 

Another possible reason for the disagreement 
with the soft pion relation is the large uncertainty 
in the matching constants. Since the vector and 
axial vector heavy-light currents are involved on 
the two sides of the equality, perturbative errors 
in the matching between the continuum and lat- 
tice operators could be important. Naively this 
error is 0(a\^, and hence should be small at /3 ~ 
6.0. Nonetheless, the large one-loop correction in 
the renormalization constant suggests that 

there could be large higher order corrections. 

In order to see how such higher order effect con- 
tributes, we computed the ratio of the renormal- 
ization constants /Zy^ nonpcrturbatively in 



the static limit, using the chiral Ward-Takahashi 
identity 



ZaZ^"^ j d^yiid^A^ - 2mqP)iy)Vo"'^ix)0) 
= -Z^^{A^\x)0), 



(1) 



where and P denote light-light axial-current 
and pseudoscalar density, with Za the renor- 
malization factor for A^, and Vi^^ and A^^ 
are the heavy-light (static-light in this particu- 
lar case) currents. We performed simulations on 
a 12'^ X 32 lattice at /3 = 6.0 following the meth- 
ods of Maiani-Martinelli ||^ . For the operator O 
we took a heavy-light meson interpolation oper- 
ator with wall source. The clover coefficient Csw 
for the light quark was chosen to be the nonper- 
turbative value from |Q 

Combining our result for Z^^ / [Zy^ Za) from 
(1) with the nonpcrturbative value of for Za, 
we obtained Zf^/Z^^ = 0.72(2) to be compared 
with the perturbative result 0.87(3). We find that 
the nonpcrturbative value for Z^^ /Zy^ is about 
20% smaller than the corresponding one-loop re- 
sult. While this explains part of the discrepancy 
between /b/Ztt and /"^((Zmax): the reduction is 
not sufficient to remove the disagreement seen in 
Fig. |. 

In our study of the renormalization constant, 
the light-light and static-light current did not in- 
clude corrections from higher dimensional opera- 
tors. Since these corrections are known to give a 
large contribution in the calculation of 0|, it 
would be important to perform a study of renor- 
malization constant with improved curents. 

4. Results for f+{q^) 

We next study the and 1 /Mb dependence of 
the form factor f^{q^). For the large q^ region, 
the form factor f^{q^) should be well approxi- 
mated by the pole dominance model with the B* 
pole, since B* is almost degenerate with B in the 
heavy quark limit and the pole is very close to 
9max- The pole model predicts 

l//+(g2) = -c,{q^ - Ml, ) + 0{{q^ - Ml,fU2) 



where the coefficient ci is written in terms of the 
dimensionless B*B7t coupling g and the B* me- 
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Figure 2. Pole fit of l//+(q^). Lattice results 
at five different heavy quark mass are shown at 
P=5.7. The B* pole is almost degenerated with 
B, and located at the origin in the plot. gMq is 
the bare heavy quark mass. 
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son decay constant /s» as ci — fn/ifB-'Mg, 
Figure || shows our data for 1/ f~^{q'^) near q^^-,^^ 
for K = 0.1369 at /3 = 5.7. We find that the pole 
fit indeed explains the data quite well for each 
value of the heavy quark mass. 

The heavy quark scaling is also predicted 
within the pole dominance model. Using fs' ^ 

1/2 3/2 

and g ~ constant, we obtain ci ~ Mg 
The slope obtained with the fit (^) is plotted 
against I /Mb in Fig. ^together with a curve rep- 
resenting We confirm that the heavy 
quark scaling is nicely satisfied. Furthermore, 
from this fit, we obtain (7=0.33(4), which is con- 
sistent with the value extracted from D* Dtt 
Isl g = 0.27(6) and with the recent lattice study 
i 5 = 0.42(8). 

5. Conclusions 

We consider that the problem of the soft pion 
relation should be resolved before a prediction of 
the B ttIV form factors from lattice QCD can 
be made. We have investigated the possibility 
that the perturbative matching contains a large 
systematic error, and found that a nonperturba- 
tive value shows a non-negligible difference from 
the one-loop result. 

Checking the dependence and the 1/Mb 
scaling of f^{q^ ) is a necessary step toward a pre- 
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Figure 3. 1/Mb scaling of the pole fit coefficient. 



diction of partial decay rate. We have found that 
the pole dominance model provides an excellent 
way to fit the observed shape of f^{q^) in the 
large q^ region. The 1/M scaling is also consis- 
tent with the pole model prediction. 
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